Neutrino Opacities at High Density and the 
Protoneutron Star Evolution 



"qo". S Reddy f [], J Pons*, M Prakash*, and J M Lattimer f 

ON 

^Department of Physics & Astronomy 
SUNY at Stony Brook, Stony Brook, NY 11794-3800, USA 

ry^ . *Departament d'Astronomia, Universitat de Valencia 

! E-46100 Burjassot, Valencia, Spain 

(N 

Abstract. The early evolution of a protoneutron star depends both on the 
CS| ' equation of state and neutrino interactions at high density. We identify the 

important sources of neutrino opacity and through model calculations show 
. that in-medium effects play an important role in determining the neutrino 

mean free paths. The effects due to Pauli-blocking and many-body correla- 
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tions due to strong interactions reduce the neutrino cross sections by large 
factors compared to the case in which these effects are ignored. We discuss 



these results in the context of neutrino transport in a protoneutron star. 

6 

H 



1. Introduction 



To date, calculations of neutrino opacities in dense matter have received relatively little 
attention jl], 0, [3], || compared to other physical inputs such as the equation of state 
(EOS). The neutrino cross sections and the EOS are intimately related. This relation- 
ship is most transparent in the long- wavelength or static limit, in which the response of 
a system to a weak external probe is completely determined by the ground state thermo- 
dynamics (EOS). Thus, in this limit, neutrino opacities consisitent with the EOS can be 
calculated [|TJ. However, when the energy and momentum transfered by the neutrinos are 
large, full consistency is often difficult to acheive in practice. Despite this, many salient 
features associated with an underlying dense matter model may be incorporated in the 
calculation of the neutrino opacities. In §2, we describe how this is accomplished for a 
non-relativistic potential model. The effects of nucleon-nucleon correlations on the neu- 
trino mean free paths are calculated using the random phase approximation (RPA) in §3, 
where we show that the magnitudes of these many-body effects on both the scattering and 
absorption reactions are large. The total scattering cross section in a multi-component 
system including the effects of correlation due to both strong and electromagnetic in- 
teractions are presented in §4 using a relativistic formalism. The implications of these 
results for neutrino transport in a protoneutron star are in §5. 

1 E-mail: reddy@nuclear.physics.sunysb.edu. 
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2. Neutrino Cross-Sections 



In the non-relativistic limit for the baryons, the cross-section per unit volume for the 
two body reaction v + B 2 — > I + B A is given by M 



G2p (V 2 + 3A 2 ) 



V 4vr 2 
The particle labels are: 1 
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incoming neutrino, 2:=incoming baryon, 3:=outgoing lepton 
(electron or neutrino), and 4:=outgoing baryon. The energy transfer to the baryons is 
denoted by q = Ei — E 3 and the momentum transfer q = \ki — k 2 \. Eq. (1) describes both 
the charged and neutral current reactions with the appropriate substitution of vector and 
axial vector coupling constants, V and A, respectively M. The factor (1 — / 3 ) accounts 
for Pauli blocking of the final state lepton. The response of the system is characterized 
by the function S(qo,q), often called the dynamic form factor. In §3, we consider the 
modifications required in the presence of spin and isospin dependent forces. 

At the mean field level, the function S(qo, q) may be evaluated exactly if the single 
particle dispersion relation is known J4|. We illustrate this using a potential model. 
In this case, if we retain only a quadratic momentum dependence, the single particle 
spectrum closely resembles that of a free gas and is given by 



Eiipi 



2M* 



n,p. 



(2) 



The single particle potentials Z7j and the effective masses M* are density dependent. 
Because the functional dependence of the spectra on the momenta is similar to that of 
the noninteracting case, it is possible to obtain an analytic expression for the dynamic 



form factor S(q , q) 



where 



Explicitly, 

S(qo,q) 



\l:;.\l:r 

nq 



exp( 
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l+exp((e ± -/i 2 + Z7 2 )/T) 
1 + exp((e± + g - + U 2 )/T) 
2q 2 



2M*x 2 




with x = 1 - (M*/M*) and c = q + U 2 - U 4 - (q 2 /2M*). The factor U 2 - U 4 is the 
potential energy gained in converting a particle of species "2" to a particle of species "4" . 
In neutral current reactions, the initial and final state particles are the same; hence, the 
strong interaction corrections are due only to M 2 . For the charged current reactions, 
modifications due to interactions are twofold. First, the difference in the neutron and 
proton single particle potentials appears in the response function and also in (i — fi n — fi p . 
Second, the response depends upon the nucleon effective masses. 

Results obtained using Eq. (3) are shown in Fig. 1, where electron neutrino absorp- 
tion (top panels) and scattering (bottom panels) mean free paths are shown in matter 
containing nucleons and leptons. Various limiting forms are also shown for compari- 
son. The charged current mean free paths are shown as thick solid lines in the bottom 
panel. These results demonstrate that effects due to kinematics, Pauli blocking, mass, 
and energy shifts are quantitatively important Q. 
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Figure 1. Absorption (top panels) and scattering (bottom panels) mean free paths in /3-equilibrated 
stellar matter. 



3. Effects of Correlations 



The random phase approximation (RPA) is particularly suited to investigate the role of 
particle-hole interactions on the collective response of matter. In this approximation, 
ring diagrams are summed to all orders. In terms of this more general response, the 
differential cross-sections are given by 



d 3 



a 



^(l - f 3 (E 3 )) R(q ,q) 
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V d 2 Q dE 3 

where H(qo, q) describes the system's response. For the neutral current reactions 
R#(qo,q) = c 2 v (1 + cos6 l )S , oo(go,g) + 4 (3 - cos 6) S 10 (q ,q) 



(4) 



(5) 



5*00 and S w are the density-density and spin-density response functions. Similarly, the 
charged current response may be written in terms of the isospin-density and the spin- 
isospin density response functions as 



Rvi(go,g)= 9 V (l + cos9)S i(q ,q) + g\ (3 - cos9)S u (q ,q) 



(6) 



3.1. Neutral Currents 

Neutrino scattering on neutrons is the dominant source of scattering opacity. We there- 
fore begin by considering the response in pure neutron matter. For a one component 
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Figure 2. Neutrino scattering cross sections versus qo/q in pure neutron matter for q—10 MeV and a 
neutrino energy of E v = 50 MeV. 



system, the RPA response functions are given by 



Sij(qo,q) 



Im n (g , q) 



-I j 



1 - exp(g /T)_ 
[1 - VijRe U°(q , q)f + [V^Im U°(q , q)} 2 . 



(7) 



The dielectric screening function is the modification introduced by the RPA to Eq. (1). 
The zeros of correspond to collective excitations, such as zero-sound and spin zero- 
sound. The potential Vy, which measures the strength of the particle- hole interaction 
in the medium, is a function of density, temperature, q , and q. Both the energy and 
momentum transferred in the particle-hole channel are of order T. Since the particle-hole 
interaction is short ranged (~l/meson mass), explicitly density dependent interactions 
play a major role in determining the magnitude of Vij. Momentum dependent interac- 
tions introduce additional structure to Vy, which in turn alters the structure of Sy. We 
will consider these latter effects separately. 

For situations in which q/k F . « 1 and T//Xj << 1, the quasi-particle interaction 
may be obtained using Fermi-liquid theory. Here, the quasi-particle interaction is given 
by = (5 2 E/5rii5nj), and is usually expressed in terms of the Fermi-liquid parameters. 
For pure neutron matter, the force in the spin independent channel is given by F , 
while the force in the spin dependent channel is given by G . Thus, Voo = F /N and 
Vio = G /N , where iVo is the density of states at the Fermi surface 0. We employ the 
results of Backmann et al [^J who have calculated F and G in pure neutron matter 
for the densities of interest here. 

In Fig. 2, the differential cross sections for v + n — > u + n are shown. The effects due 
to correlations lead to significant reductions at small q^, since the interaction is repulsive 
and the excitation of collective modes enhances the response at large go- The well defined 
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Figure 3. Charged current reaction cross sections in symmetric nuclear matter. 



zero-sound and spin zero-sound seen at T = are damped at finite temperatures. The 
contribution from the low q Q region dominates the total cross-section due to final state 
blocking. The presence of a repulsive particle-hole force acts to reduce the neutrino cross 
sections for neutrino energies of order T. 

The composition of charge neutral, /3-equilibrated stellar matter depends on the 
EOS of dense matter, which contains an admixture of protons and electrons. The cor- 
relations due to interactions between the different particle species and electromagnetic 
correlations also play an important role. In §4, we discuss a relativistic framework to 
describe the RPA response of a mixture of neutrons, protons, and electrons. 

3.2. Charged Currents 

The charged current reaction is kinematically different from the neutral current reaction, 
since the energy and momentums transfers are not limited only by the matter tempera- 
ture. The energy transfer is typically of order fi = \i n — \i p . The charged current probes 
smaller distances; hence, the go and q dependencies of the particle-hole force are likely to 
play an important role. In symmetric nuclear matter jl = 0, and therefore, the situation 
is very similar to the neutral current case. Since the collective response of laboratory 
nuclei is well established in nuclear physics, we begin by considering the response of sym- 
metric nuclear matter to a charged current weak probe. Here the strength parameters 
Voi for the spin independent iso-spin channel and the V\\ for the spin dependent iso-spin 
channel are required. In nuclear matter, the particle-hole force (retaining only the I = 
terms) is given by || 

F(h, k 2 ) = N^[F + GVi<x 2 + nr 2 (^ + G' Q {a x a 2 )\ . (8) 

For the charged current reaction, isospin and charge are transferred along the particle- 
hole channel; hence, only the last two terms in Eq. (8) contribute. The potentials required 
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to calculate the RPA response in Eq. (7) are given by Vio = 2Fq/N and Vn = 2G' /N , 
where the factor two arises due to isospin considerations. The Fermi-liquid parameters 
may be calculated from an underlying dense matter model, as for example, a Skyrme 
model, which successfully describes the low-lying excitations of nuclei @. In Fig. 3, the 
response of symmetric nuclear matter to the charged current probe is shown. Results are 
for two different parameterizations of the Skyrme force, SGII and SkM*. The excitation 
of giant-dipole and Gamow- Teller resonances shifts the strength to large qo, and as in 
the case of scattering, the low go response is significantly suppressed. 
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Figure 4. Neutrino absorption cross sections in asymmetric matter. 

In asymmetric matter, the energy and momentum transfers are large. Therefore, 
the momentum dependence of the particle-hole force becomes important. At high mo- 
mentum transfer, the conservation of the vector current implies that the response func- 
tion Sio is not strongly modified. Further, since its contribution to the differential cross 
section is roughly three times smaller than the spin-dependent response function Sn, 
we focus on the Gamow- Teller part and assume that the Fermi matrix elements are not 
screened. The iso-vector interaction in the longitudinal channel arises due to 7r exchange, 
and in the transverse channel due to p meson exchange. In addition, to account for the 
large repulsion observed, Migdal || introduced screening in this channel through the 
parameter g' . This form for the particle-hole interaction has been successful in describ- 
ing a variety of nuclear phenomena. The longitudinal and transverse potentials in the 
17 + P + 9' model are given by H 



v L (qo,q) = 

V T (q , q) = 
where F n = (A 2 - m 2 r )/(A 2 - q 2 ) and F f 




^(q)+g' 



2 F 2 p (q)+9' 



(9) 



(A 2 - m 2 ) /(A 2 - q 2 ) are the ttNN and pNN 
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form factors. Numerical values used are C p — 2, g' — 0.7, A = 1.4 GeV, and A p = 2 GeV. 
The RPA response function then takes the form 



Su{qo,q) 



exp(g /T) 



/ m n» to , ,)(- + - 



1 



e L , T = [1 -2V L>T ReU°(q ,q)} 2 + [2V LiT Im U°(q , q)f 



(10) 



The cross sections in charge neutral stellar matter with a fixed lepton fraction Yl = 0.4 
are shown in Fig. 4. The RPA screening reduces the cross sections by about a factor of 
two. The inclusion of virtual A-hole excitations reduces the cross sections further. Unlike 
in the case of ordinary f3 decay, where the A-hole contribution is important, the role of 
these excitations is marginal here due to the large energy and momentum transfers. 



4. Relativistic Treatment 




Figure 5. Neutrino scattering cross sections in a relativistic approach. 

At high density the baryons become increasingly relativistic and thus a relativis- 
tic description may be more appropriate. Such a description also allows us to treat the 
baryons and the electrons on an equal footing. The basic formalism and some illustrative 
results may be found in |3|, §. Here we present results for a system of neutrons, protons, 
and electrons at finite temperature. The ground state is described by a field-theoretical 
model at the mean field level, in which isoscalar a, uo and isovector p meson exchanges 
modify the in-medium baryon propagators. We calculate the Hartree response by ac- 
counting for these modifications. Correlations are incorporated through the relativistic 
random phase approximation (RRPA), where the particle- hole interaction is mediated 
by the a and u mesons between the baryons, and electromagnetic interactions between 
protons and electrons are mediated by photons. Fig. 5 shows the results for different 



7 



temperatures at a baryon density of 0.32 fm -3 . Compared to the non-relativistic models, 
the screening effects are small, the dominant suppression arising mainly due to density 
dependent nucleon effective masses. 

5. Discussion 

We have highlighted the influence of correlations and collective phenomena on the neu- 
trino opacities in dense matter. Our findings here indicate that the neutrino cross sections 
are significantly reduced, and the average energy transfer in neutrino-nucleon interac- 
tions is increased due to the presence correlations in the medium. Several improvements 
are necessary before we can assess the influence of these results on the macrophyscial 
evolution of a protoneutron star. Among the most important of these are calculations 
that provide (1) the dynamic form factor, (2) the particle-hole and particle-particle in- 
teractions, (3) the renormalization of the axial charge, and, (4) the means to assess the 
role of multi-pair excitations, in charge neutral, beta-equilibrated dense matter at finite 
temperature. While investigations along these directions are in progress (see also U), 
some general trends may be anticipated. The many-body effects studied here, including 
the improvements listed above, suggest considerable reductions in the opacities compared 
to the free gas estimates often employed in many applications. In the particular instance 
of the early evolution of a protoneutron star, the suggested modifications imply shorter 
time scales over which the deleptonization and cooling occur. Since mean free paths A 
are much less than the stellar radius R, evolution is via diffusion. To order of magni- 
tude, therefore, the timescale r oc 1/A. However, there are important feedbacks between 
these evolutionary consequences and the underlying EOS, in particular the specific heat 
of multicomponent matter, which have to be studied before firm conclusions may be 
drawn. 
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